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Abstract 

We describe a natural strategy to enumerate compact hyperbolic 3- 
manifolds with geodesic boundary in increasing order of complexity. 
We show that the same strategy can be employed to analyze simul- 
taneously compact manifolds and finite- volume manifolds having toric 
cusps. In opposition to this we show that, if one allows annular cusps, 
the number of manifolds grows very rapidly, and that our strategy can- 
not be employed to obtain a complete list. We also carefully describe 
how to compute the volume of our manifolds, discussing formulae for 
the volume of a tetrahedron with generic dihedral angles in hyperbolic 
space. 

MSC (2000): 57M50. 

According to Thurston's geometrization program, the theory of hyperbolic 
manifolds plays a central role in 3-dimensional topology. Hyperbolic man- 
ifolds with geodesic boundary, the first example of which was given by 
Thurston himself in jU] (and later generalized in (101), are an important 
portion of this theory. On the other hand, the algorithmic and computer 
approach to 3-manifolds has been acquiring an increasing popularity in 
recent years. For cusped hyperbolic manifolds this approach, which was 
worked out in [3J 023 02] and several other papers, again depends on 
ideas of Thurston, namely on the use of moduli for hyperbolic ideal tetra- 
hedra and equations to ensure consistency of the structures. For closed 
manifolds the basics of the computer approach were set by Matveev in |34| 
(see also ^JJEHl)) and several experimental results were later obtained by 
himself and other authors (see |3H I32| ). In the present paper we describe 
the general setting of the algorithmic approach to hyperbolic 3-manifolds 

*Both authors were supported by the INTAS Project "CalcoMet-GT" 03-51-3663 
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with geodesic boundary, concentrating in particular on their enumeration in 
order of increasing complexity, and on the computation of their volume. 

1 Hyperbolic structures 

In this section we review the general theory of hyperbolic 3-manifolds with 
geodesic boundary, stating the main results we will need in the sequel. 

Local structure, cusps, compact ificat ions In the rest of this paper 
we will call hyperbolic an orientable finite-volume complete Riemannian 3- 
manifold with non-empty boundary, locally isometric to an open subset of a 
closed half-space of H 3 . We will always denote such a manifold by Y. Note 
that dY is totally geodesic. Doubling Y along its boundary and using the 
description of the ends of the double [2], one can show that Y consists of 
a compact portion together with some "toric and annular cusps." A toric 
cusp is here a space of the form T x [0, oo), attached to the rest of Y along 
T x {0}, and an annular cusp is defined analogously. Note that toric cusps 
are disjoint from dY, while an annular cusp gives two punctures in dY. In 
particular, dY is compact if and only if Y has no annular cusps. 

Given Y as above, we can naturally get a compact manifold Y by adding 
a torus T x {oo} for each toric cusp T x [0, oo), and an annulus A x {oo} 
for each annular cusp A x [0, oo). However, it turns out that when there 
are annular cusps another compactification Y of Y is more suited to the 
geometric situation. We define Y as a quotient of Y, where in each annular 
cusp Ax[0, oo] with A = S 1 x [0, 1], for all t € [0, 1] we collapse S 1 x {t} x {oo} 
to one point {*} x {t} x {oo}. Note that Y is obtained from Y by removing 
some toric boundary components and drilling some properly embedded arcs. 

Rigidity and Kojima decomposition A key result for computational 
purposes is the following: 

Theorem 1.1. Any two homeomorphic hyperbolic manifolds are isometric. 

This result is commonly referred to as rigidity theorem, and a proof for 
the geodesic boundary case was spelled out in [17) . 

Another very important fact is that the hyperbolic structure determines 
certain combinatorial data which can be employed to efficiently test two 
manifolds for homeomorphism. This result is analogous to the Epstein- 
Penner decomposition of cusped hyperbolic manifolds without boundary 12 , 
and it was proved by Kojima |26| I27j . Its statement involves the notion of 
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truncated polyhedron, that we now give. Consider the projective model of 
hyperbolic 3-space, viewed as the open unit ball B 3 in the Euclidean 3-space 
E 3 . Let us call finite the points of B 3 , ideal those of 9B 3 , and ultra-ideal the 
other points of E 3 . Consider a convex polyhedron P in E 3 with ideal and/or 
ultra-ideal vertices, and all edges meeting the closure of B 3 . Dual to each 
ultra-ideal vertex of P there is an open hyperbolic half-space, and we define 
P to be P n B 3 minus these half-spaces. Any P arising like this from some 
P will be called a truncated polyhedron. Note that P has internal faces, 
those coming from faces of P, and truncation faces. Moreover internal and 
truncation faces lie at right angles to each other. 

Theorem 1.2. Any hyperbolic manifold admits a canonical decomposition 
as a gluing of truncated polyhedra along the internal faces. 

In the sequel we will need to refer to the geometric argument underlying 
this result, so we briefly sketch it here. Regard E 3 as the hyperplane at 
height 1 in Minkowsky 4-space M 3,1 . Suppose first that Y has no toric 
cusps, and identify the universal cover of Y to an intersection of half-spaces 
of B 3 . Dual to each such subspace there is a point having norm 1 in M 3,1 , 
and the decomposition of Y is obtained by taking the faces of the convex hull 
of all these points, projecting first to E 3 , truncating, and then projecting to 
Y. When there are toric cusps one must also take suitably small Margulis 
neighbourhoods of these cusps, lift them to horoballs in B 3 , consider the 
duals to these horoballs on the light-cone of M 3 ' 1 , include these duals in the 
convex hull, and suitably subdivide some of the resulting faces. We address 
the reader to |17j for all the details. 

We are now in a position to explain why we have introduced the com- 
pactification Y: 

Proposition 1.3. Let Y be hyperbolic. For each truncated polyhedron P 
in the Kojima decomposition of Y , consider the corresponding Euclidean 
polyhedron P. Glue these P's along the same maps as in the Kojima de- 
composition, and remove open stars of the vertices. The resulting space is 
then homeomorphic to Y . 

Topological obstructions to hyperbolicity Since dealing with open 
manifolds is impossible by computer, the idea to enumerate the hyperbolic 
Y f s is to enumerate the corresponding compact Y's or Y's. The next two 
results show that the presence of annular cusps makes a dramatic difference. 
The former can be found in ^7], the latter follows from the results discussed 
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in Section |HJ (Recall that we are calling 'hyperbolic' a manifold with non- 
empty geodesic boundary). 

Theorem 1.4. Let M be a compact orientable 3-manifold with boundary, 
and let Y be M minus the toric components ofdM. The following conditions 
are pairwise equivalent: 

• Y is hyperbolic; 

• Y is hyperbolic, it has no annular cusps, and Y = M ; 

• M is irreducible, boundary-irreducible, acylindrical and atoroidal, and 
X (M) < 0. 

Proposition 1.5. • For any g 2 there exists a hyperbolic Y such that 
Y is the handlebody of genus g; 

• For any compact 3-manifold M there exists a hyperbolic Y such that 

2 Complexity of 3-manifolds 
and enumeration strategy 

In this section we recall the basics of the theory of simple and special spines, 
and the related theory of ideal triangulations, describing how it can be 
employed to (partially) enumerate the class of 3-manifolds we are interested 
in. Proofs of all results on spines and complexity can be found in |35j . 

Simple spines and complexity Throughout the present section we will 
employ the PL category for 3-manifolds and use the customary notions of 
PL topology, see [321- A simple polyhedron is a compact polyhedron P such 
that the link of each point of P can be embedded in the space given by a 
circle with three radii. In particular, P has dimension at most 2. Finite 
graphs and closed surfaces are examples of simple polyhedra. A point of a 
simple polyhedron is called a vertex if its link is precisely given by a circle 
with three radii. A regular neighbourhood of a vertex is shown in Fig. ^ 
(3). From the figure one sees that the vertices are isolated, whence finite in 
number. Graphs and surfaces do not contain vertices. 

If M is a compact 3-manifold with non-empty boundary, we call spine 
of M a subpolyhedron P of M such that M\P is an open collar of dM. We 
call complexity of M, and denote by c(M), the minimal number of vertices 
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Figure 1: Local aspect of an almost-special polyhedron. 



of a simple spine of M. We say that a spine of M is minimal if it has c(M) 
vertices and it does not contain any proper subpolyhedron which is also a 
spine of M. 

Special spines and ideal triangulations We now introduce a more 
restrictive type of spine which turns out to have a very clear geometric 
counterpart. A simple polyhedron P is called almost-special if the link of 
each point of P is given by a circle with either zero, or two, or three radii. 
The local aspects of P are correspondingly shown in Fig. ^ The points of 
type (2) or (3) are called singular, and the set of singular points of P is 
denoted by S(P). We will say that P is special if it is almost-special, S(P) 
contains no circle component, and P \ S(P) consists of open 2-discs. 

We now call ideal triangulation of a compact 3-manifold M with non- 
empty boundary a realization of the interior of M as follows. We take a finite 
number of tetrahedra, we glue together in pairs the faces of these tetrahedra 
along simplicial maps, and we remove the vertices. Equivalently, an ideal 
triangulation of M is a realization of M as a gluing of truncated tetrahedra. 
The relation between spines and triangulations is given by the following: 

Proposition 2.1. The set of ideal triangulations of a 3-manifold M corre- 
sponds bijectively to the set of special spines of M . The polyhedron corre- 
sponding to a triangulation is the 2-skeleton of the dual cellularization, as 
shown in Fig. 

Manifolds having special minimal spines Special spines have two 
main advantages if compared to merely simple ones. First of all, a special 
spine determines the manifold it is a spine of [3], which is false for sim- 
ple spines. Second, no efficient method for listing simple spines is known, 
whereas enumerating special spines in increasing order of complexity is very 



5 



Figure 2: Duality between triangulations and special polyhedra. 



easy (at least theoretically, see [3]). For these reasons, the next result is 
crucial for computational purposes: 

Proposition 2.2. Let M be a compact 3-manifold with non-empty bound- 
ary. The following conditions are pairwise equivalent: 

• M is irreducible, d -irreducible, acylindrical, and not the 3-disc; 

• M has some special minimal spine; 

• All minimal spines of M are special. 

This result is not quite stated in this form in or in any of Matveev's 
papers, but it easily follows from the proof of |35| Theorem 2.2.4]. 

Naif enumeration strategy Cusped hyperbolic manifolds without boun- 
dary were studied and enumerated in [I], which explains why we have de- 
cided to restrict to manifolds with non-empty boundary. In addition, we 
forbid here the presence of annular cusps, because, according to Proposi- 
tions ED an d E21 the theory of spines, triangulations, and complexity does 
not appear to be well-suited for the investigation of such manifolds. See 
Section H3 

Let us then define TL n as the set of all hyperbolic manifolds having com- 
plexity n and non-empty compact geodesic boundary. We also define TL n as 
the set of all orientable, compact, irreducible, <9-ir reducible, and acylindri- 
cal manifolds with negative \- Theorem 11.41 implies that if Y £ TL n then 
Y € 7i n . Conversely, if M G 7i n , then M minus the toric components of dM 
belongs to ri n if and only if M is atoroidal. We can then view 7i n as the set 
of candidate hyperbolic manifolds with complexity n and without annular 
cusps. 
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Note now that Proposition 12.21 applies to the elements of 7i n . There- 
fore the following theoretical steps lead to the exact determination of 7i n , 
assuming rt m is known inductively for m < n: 

1. Produce the list of all special spines with n vertices and negative x; 

2. Remove from the list the spines P whose associated manifold M(P) is 
not hyperbolic; 

3. Remove from the list the spines P such that M(P) belongs to some 
TC m for m < n; 

4. For P varying in the list, compare the manifolds M(P) for equality, 
discarding duplicates. 

Step (1) does not present any theoretical difficulty, but its practical 
implementation is quite demanding if no computational shortcuts are em- 
ployed. We will discuss these shortcuts in the next paragraph. And we will 
explain how to carry out the other steps in the next section. 

Pseudo-minimal spines Taking into account Propositions 12.11 and 12.21 

the reader may wonder why we have employed special spines rather than 
triangulations to list the elements of Ti n . The first remark is that one can 
always restrict to minimal triangulations without losing any potentially in- 
teresting manifold. However, it turns out that dual to a triangulation which 
is minimal among triangulations of the same manifold, there is often a spe- 
cial spine which is not minimal among simple spines of the same manifold. 
Of course this can only happen if the corresponding manifold violates one 
of the topological constraints of Proposition 12.21 but we know from Theo- 
rem ^31 that in this case the manifold is not hyperbolic, so we can discard 
it. In other words, non-minimality is a much more flexible notion for spines 
than for triangulations, so using spines we can substantially reduce the list 
of manifolds that we will later need to investigate. 

Of course it is impossible to check in a direct fashion whether a special 
spine is minimal, but there are many criteria for non-minimality, which can 
be used as tests to discard spines which will certainly not bring any relevant 
manifold. The tests used in JB] are based on the moves shown in Fig. |31 
which are easily seen to transform a spine of a manifold into another spine of 
the same manifold. To be precise, let us say that a spine is pseudo-minimal if 
it cannot be transformed into a spine with fewer vertices by a combination of 
the moves shown in Fig. 01 The key point mentioned above on the flexibility 
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Figure 3: Moves on simple spines. 



of spines is that moves (1) and (2) do not lead to special spines (in general), 
so they do not have counterparts at the level of triangulations. The first 
step of the enumeration strategy is then replaced by the following: 

• Produce the list of all pseudo-minimal special spines with n vertices 
and negative x- 

We also mention that another trick very important for computational 
purposes is to construct the candidate spines portion after portion, following 
the branches of a tree, and to apply the non- minimality tests arising from 
the moves of Fig. 01 also to partially constructed spines, thus "cutting the 
dead branches" at their bases. 

3 Hyperbolicity equations and tilts 

In this section we discuss how an ideal triangulation can be employed to 
construct a hyperbolic structure on a given manifold and to recognize the 
canonical Kojima decomposition of that manifold. This allows to carry out 
steps (2)- (4) of the enumeration strategy for 7i n explained in the previous 
section. We actually include in the discussion also manifolds with annular 
cusps, because the methods via ideal triangulations to construct and recog- 
nize the hyperbolic structure apply to these manifolds too. It is only step 
(1) of the enumeration strategy (the listing of special spines) that is not 
suited to manifolds with annular cusps, and the reason why these manifold 
are ruled out from 7i n . 

We first treat the compact case and then sketch the variations needed 
for the case where there are also some cusps. For all details and proofs (and 
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for some very natural terminology that we use here without giving actual 
definitions) we address the reader to |17j . 

Compact case: moduli and equations The basic idea for constructing 
a hyperbolic structure via an ideal triangulation is to realize the tetrahedra 
of the triangulation as truncated tetrahedra in H 3 and then require that 
the structures match when the tetrahedra are glued together. The following 
facts show that one can use the dihedral angles as moduli to parameterize 
the realizations of a tetrahedron and to check consistency: 

• A hyperbolic structure on a combinatorial truncated tetrahedron is 
determined by the 6-tuple of dihedral angles along the internal edges; 

• The only restriction on this 6-tuple of positive reals is that the angles 
of each of the four truncation triangles should sum up to less than tt; 

• The lengths of the internal edges can be computed as explicit functions 
of the dihedral angles; 

• A choice of hyperbolic structures on the tetrahedra of an ideal trian- 
gulation of a manifold M gives rise to a hyperbolic structure on M 
if and only if all matching edges have the same length and the total 
dihedral angle around each edge of M is 2tt. 

Given an ideal triangulation consisting of n tetrahedra one then has the 
hyperbolicity equations: a system of 6n equations with unknown varying in 
an open set of R 6n which, by rigidity (Theorem II. 1JI . admits one solution at 
most. 

Canonical decomposition and tilts Once a hyperbolic structure has 
been constructed on a manifold Y using an ideal triangulation 7", one natural 
issue is to decide if T is the canonical decomposition of Y and, if not, to 
promote T to become canonical. These matters are faced using the tilt 
formula (Ml EDI ESI, that we now describe. 

Recall first that the Kojima decomposition of Y is constructed by pro- 
jecting first to H 3 and then to Y the faces of the convex hull of the set 
V C M 3 ' 1 of the duals to the boundary components of the universal cover 
of Y. If o" is a ci-simplex in T, each end of a lifting of a to H 3 determines 
a point of V . Now let two tetrahedra and A2 share a 2-face F, and let 
Ai, A 2 and F be liftings of A x , A 2 and F to H 3 such that Ai n A 2 = F. 
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Let F be the 2-subspace of M • that contains the three points of V deter- 

~ rp\ 

mined by F. For i = 1, 2 let A i be the half-3-subspace bounded by F 

and containing the fourth point of V determined by Aj. Then one can show 

that T is the canonical Kojima decomposition of Y if and only if, whatever 

F, Ai,A2, the following conditions are met: 

(a) the convex hull of A^ and A% does not contain the origin of M 3,1 ; 

(b) A^ and A 2 F) lie on distinct 3-subspaces of M 3>1 . 

In addition, if condition (a) is met for all triples .F, Ai,A 2 , the canonical 

decomposition is obtained by merging together the tetrahedra along which 

condition (b) is not met. 

The tilt formula defines a real number t(A, F) describing the "slope" 
/ p\ 

of A . More precisely, one can translate conditions (a) and (b) into the 
inequalities t(Ai, F)+t(A 2 , F) < and i(Ai, F)+i(A 2 , F) / respectively. 
Since we can compute tilts explicitly in terms of dihedral angles, this gives a 
very efficient criterion to determine whether T is canonical or a subdivision 
of the canonical decomposition. Even more, it suggests where to change T 
in order to make it more likely to be canonical, namely along 2-faces where 
the total tilt is positive. This is achieved by 2-to-3 moves along the offending 
faces, as discussed in [T7] , 

The non-compact case When one is willing to accept both compact 
geodesic boundary and cusps, the same strategy for constructing the struc- 
ture and finding the canonical decomposition applies, but many subtleties 
and variations have to be taken into account. Let us quickly mention which. 

Moduli. As suggested by Proposition ll.31 to construct a hyperbolic struc- 
ture on a manifold Y one must take an ideal triangulation T of the com- 
pactification Y, in which each toric cusp is completed with a torus and each 
annular cusp is completed with a segment. Moreover one must suppose 
that the segments arising from the annular cusps of Y are edges of T. To 
parameterize tetrahedra one should then proceed as follows: 

• Assign dihedral angle to each edge corresponding to an annular cusp 
of Y, which geometrically means that the edge, before truncation, is 
tangent to the boundary of B 3 ; 

• If three edges meet at a vertex asymptotic to a toric cusp of Y, assign 
them dihedral angles summing up to 7r, which geometrically means 
that the vertex is an ideal one. 
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Switching viewpoint, if one starts from an ideal triangulation T of a 
manifold M, one must arbitrarily choose a family a of edges of T and assign 
dihedral angles to the edges in a, and dihedral angles summing up to ir 
to triples of edges asymptotic to toric components of 8M. If the consis- 
tency and completeness equations are satisfied (see below), this leads to a 
hyperbolic structure on M minus a and the toric components of dM. 

Equations. If an internal edge with non-zero dihedral angle ends in a 
cusp then its length is infinity, so some of the length equations must be 
dismissed when there are cusps. There are no consistency issues connected 
with half-infinite edges but, when an edge is infinite at both ends, one must 
make sure that the gluings around the edge do not induce a sliding along the 
edge, which translates into the condition that the similarity moduli 2 of the 
Euclidean triangles around the edge have product 1. This ensures existence 
of the hyperbolic structure, but one still has to impose completeness of cusps. 
Just as in the case where there are cusps only, this amounts to requiring that 
the similarity tori on the boundary be Euclidean, which translates into the 
holonomy equations involving the similarity moduli. Note that there is no 
completeness issue connected to annular cusps. 

Canonical decomposition. When there are cusps, the set of points V to 
take the convex hull of consists of the norm-1 duals of the boundary com- 
ponents of the universal cover and of some points on the light-cone dual 
to Margulis neighbourhoods of the cusps. The precise discussion of how to 
choose these extra points is too complicated to be reproduced here (see |17j). 
but in practice one has that the choice of sufficiently small Margulis neigh- 
bourhoods always leads to the right result. 

4 Volume computation 

All the hyperbolic manifolds found by computer in can be decomposed 
into genuine tetrahedra, i.e. tetrahedra with positive dihedral angles (with 
ultra-ideal vertices, and possibly some ideal ones). We do not know if all 
hyperbolic manifolds admit such a genuine decomposition, but we know 
that if we subdivide the Kojima decomposition into tetrahedra we find some 
genuine tetrahedra and possibly some flat ones, which do not contribute to 
the volume. Therefore the problem of computing the volume of a hyperbolic 
manifold is completely reduced to the same problem for a single tetrahedron 
with arbitrarily assigned dihedral angles. We discuss in this section various 
formulae for this volume. In doing this we will employ a unified viewpoint, 
which includes finite, ideal, and ultra-ideal vertices. 
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Volumes of polyhedra: historical remarks The calculation of the 
volume of a polyhedron in 3-dimensional space is a very old and difficult 
problem. The first known result in this direction belongs to Tartaglia (1494) 
who found a formula for the volume of a Euclidean tetrahedron, now known 
as the Cayley-Menger formula. It was recently shown in [33] and that 
the volume of any Euclidean polyhedron is a root of an algebraic equation 
whose coefficients are polynomial functions of the lengths of the edges, the 
polynomials being determined by the combinatorial type of the polyhedron. 

In the hyperbolic and spherical spaces the situation is much more com- 
plicated. We confine here to the hyperbolic case, and we first recall that 
the volume formula for a biorthogonal tetrahedron (orthoscheme) has been 
known since Lobachevsky and Schlafli (see |3L)j and jlB] respectively). The 
volumes of the Lambert cube and of some other polyhedra were computed 
by Kellerhals |24| . Derevnin and Mednykh [S], Mednykh, Parker, and Ves- 
nin [HE], and others. The volume of a regular tetrahedron in hyperbolic 
space was investigated by Martin [33]. The volume formula for a hyperbolic 
tetrahedron with a few non-ideal vertices was found by Vinberg |52| . 

Despite these partial results, a formula for the volume of an arbitrary 
hyperbolic tetrahedron has been unknown until very recently. The general 
algorithm for obtaining such a formula was indicated by Hsiang [22] , and the 
complete solution of the problem was given within the space of a few years 
by several authors, namely Cho and Kim _6 , Murakami and Yano |39j . and 
Ushijima [SJ. In all these papers the volume of a tetrahedron is expressed 
as an analytic formula involving 16 dilogarithm or Lobachevsky functions 
whose arguments depend on the dihedral angles of the tetrahedron and on 
some additional parameter which is found as a root of some complicated 
quadratic equation with complex coefficients. 

A geometric meaning of the Murakami- Yano formula was recognized 
by Leibon [2H] from the viewpoint of the so-called Regge symmetry. An 
excellent exposition of these ideas and a complete geometric proof of the 
Murakami- Yano formula was given by Mohanty |38| . It is worth mentioning 
that the ideas of Regge symmetry and scissors congruence were also partially 
used by Cho and Kim in [£], where the first general formula was actually 
obtained. 

A remarkable phenomenon is that the volume formulae for the tetra- 
hedron become much easier for symmetric tetrahedra, i.e. for tetrahedra 
with identical dihedral angles at opposite edges. This fact was first noticed 
by Milnor [37], who expressed the volume of an ideal tetrahedron (which is 
automatically symmetric) as the sum of the values of the Lobachevsky func- 
tion on the three dihedral angles. It was later shown by Derevnin, Mednykh, 
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Figure 4: Notation for the dihedral angles of a tetrahedron. 



and Pashkevich that a rather simple formula also exists for an arbitrary 
symmetric tetrahedron with finite vertices. 

In the next paragraph we will present an elementary integral formula for 
the volume of an arbitrary tetrahedron in hyperbolic space. The formula 
involves some parameters depending on the dihedral angles, and it is very 
helpful to actually evaluate the volume by computer. The Murakami- Yano 
result can be obtained as an easy consequence of this formula. 

Volume of a tetrahedron without truncation Let T(A, B, C, D, E, F) 
denote the tetrahedron in hyperbolic 3-space with dihedral angles A, B, C, 
D, E, F as in Fig. |1J Recall ^Z] that T is determined up to isometry by 
A, B, C, D, E, F. Moreover if the sum of the triple of angles at a certain 
vertex is ir then the vertex is ideal. Similarly, if the sum is less than tt then 
the vertex is ultra-ideal, i.e. the tetrahedron is truncated at the vertex, 
which is then replaced by a triangle. 

The next result is due to Derevnin and Mednykh |1U| : 

Theorem 4.1. Suppose that all the vertices ofT = T(A, B, C, D, E, F) are 
ideal or finite. Then Vol(T) is given by 




B+D+F+z 
2 



• sm 



B+C+E+F+z 
2 



• COS 



C+D+E+z 
2 



•sinf 



dz 
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with z\ and z 2 given by 




h ' k s ' 




where 



k x = -(cos S + cos(A + D) + cos{B + E) + cos(C + F) 
+ cos(F + E + F) + cos(F + B + C) 
+ cos(^ + E + C) + cos(,4 + B + F)), 

£) 2 = sin S 1 + sin(A + D) + sin(5 + E) + sin(C + F) 
+ sin(F) + F + F) + sin(F + B + C) 
+ sin(^ + F + C) + sin(A + F + F), 

Ai3 = 2(sinAsinF + sin F sin F + sin C sin F), 



and S = A-\-B + C + D + E + F. Moreover the kj 's and Zj 's are all real 
numbers, so the integral is just an ordinary integral on an interval of the 
real line, and the function to be integrated vanishes at the Zj 's. 

Remark 4.2. There is a very transparent geometric interpretation of the 
sums of dihedral angles V\ = A+B+C, V 2 = A+E+F, V 3 = B+D+F, and 
V4 = C + D + F appearing in the numerator of the volume formula, and of 
the sums H x = A+B + D + E, H 2 = A+C + D + F, and F 3 = B + C + E + F 
appearing in the denominator. Namely, the Vj's correspond to the triples of 
edges incident to the vertices of the tetrahedron, while the Fj's correspond 
to the Hamiltonian cycles. 

Remark 4.3. The parameters z\ and z 2 appearing in the volume formula 
can be shown to be roots of the equation k\ cosz + k 2 s'mz = k%, while 
k\ = — 4det(G), where G is the Gram matrix of T. The numbers z\ and 
z 2 also have a geometric meaning, as explained in j^H]. Namely, they arise 
as parameters for the decomposition of an ideal octahedron into four ideal 
tetrahedra with an edge in common. The octahedron is canonically defined 
by the tetrahedron T, and its dihedral angles are just linear combinations 
of those of T. 

Recall now that the dilogarithm function is defined by the integral 




X 
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where x € C \ [1, oo) and log is the continuous branch of the logarithm 
function given by log £ = log |£| +i arg £ with the constraint — tt < arg £ < ix. 
Let 

l(z) = U 2 (e iz ). 

As an immediate consequence of Theorem l4. ll we have the following Murakami- 
Yano-Ushijima formula obtained in |39j and |51j : 

Corollary 4.4. Suppose that all the vertices ofT = T(A, B, C, D, E, F) are 
ideal or finite. Then 



Vol(T) = ^(u(z 1 ,T)-U(z 2 ,T) 



where 



U(z,T) = ~(l( z ) + l(A + B + D + E + z) 

+l(A + C + D + F + z) + l{B + C + E + F + z 
-Z(tt + A + B + C + z) -l(ir + A + E + F + z) 

— Z(tt + B + D + F + z)-l{ir + C + D + E + z) 

Conerning this formula we note that 

= 9f(Li 2 (e")) = 2A(z/2), 

where A(z) is the Lobachevsky function defined by the integral 

e 

A(6») = - j log|2sint|dt. 
o 

This shows, as already announced, that the volume of a tetrahedron is an 
algebraic sum of 16 Lobachevsky functions. 

Truncated tetrahedra According to |51| and |52| . the determinant of the 
Gram matrix is strictly negative also for truncated tetrahedra, so the num- 
bers fei, . . . , &4 appearing in the statement of Theorem l4,ll are still real. How- 
ever a non-trivial issue arises concerning the choice of the analytic branch of 
the arctan function for the definition of the numbers z\ and z 2 - More pre- 
cisely, one can be forced to change branch, in order to ensure the continuity 
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of Z\ and Z2 , when k\ approaches 0, which indeed can happen when there 
are some ultra-ideal vertices. For example, for 

( 7T 7T TT TT TT 7T \ 

(A, B, C, D, E,F)= — , -, , — , -, 

v ' ' ' ' ' ' V12 3 10.18 12 3 10.18/ 

we have k\ > 0, whereas for 

„ N / TT TT TT TT TT TT \ 

(A, B, C, D, E,F)= — , -, , — , -, 

v ' ' 1 V12 3 10.19 12 3 10.19/ 

we obtain k\ < 0. Therefore in this case we have to change the analytic 
branch of arctan^/^i)- The most convenient way to do so is to replace 
arctan(/c2 Ai) by tt/2 — arctan(fci/fc2), which yields a continuous variation 
of the values of z\ and Z2 as k\ passes through 0. We note that k% and k^ 
are always positive, so there is no such problem with arctan(/c4//c3). One 
can now check that, after appropriately choosing the analytic branches of 
the arctan function, the formula of Theorem 14. II still gives the correct value 
of the volume for truncated tetrahedra. 



Symmetric tetrahedra We will say that T = T(A, B, C, D, E, F) is sym- 
metric if A = D, B = E, and C = F, and in this case we will denote T just 
by T(A, B,C) for simplicity. 

We begin with the easy case of ideal tetrahedra, which are automatically 
symmetric. Moreover T(A, B, C) is ideal if and only if A + B + C = it. In 
this case the volume turns out [31! to be given simply by 

Vol(T(A B, C)) = A(A) + A(B) + A(C), 

where A is the Lobachevsky function defined above. For symmetric tetra- 
hedra T(A,B,C) with finite vertices, i.e. such that A + B + C > tt, the 
following was shown in |llj : 

Theorem 4.5. If T = T(A, B, C) has finite vertices then Vol(T) is given 
by 

tt/2 

f sin^ 1 (cos A cos t) + sin -1 (cos B cost) + sin -1 (cos C cost) — sin -1 (cost) 

2 / 1 

J sin It 

e 

where 6 6 (0, tt/2) satisfies 

1 - a 2 - b 2 - c 2 - 2abc 

tanO = — 

a/(1 - a + b + c)(l + a - b + c)(l + a + b - c)(-l + a + b + c) 
with a = cos A, b = cos B, and c = cos C. 
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Remark 4.6. The value of 9 in the previous statement has a simple ge- 
ometric interpretation in terms of the "sine rule" (see Theorem 7]): 
if Ia,Ib,Ic are the lengths of the edges of T with dihedral angles A,B,C 
respectively, then 

sinj4 sinB sinC ^ 
sinh I a sinh I b sinh Iq 

We conclude this section by noting that no simple formula is currently 
known for the volume of a symmetric tetrahedron with ultra- ideal vertices, 
but it is reasonable to expect that such a formula should exists and have 
connections with an ultra-ideal version of the sine rule. 



5 Manifolds without annular cusps 

In this section we recall the definition and main properties of a certain 
class of manifolds M. g ^ C H-g+k studied in and we describe the 

experimental results of JH] on 7i n for n ^ 4. Recall that TL n contains 
the hyperbolic manifolds of complexity n with non-empty compact geodesic 
boundary. The next section will be devoted to manifolds with non-compact 
boundary, i.e. with annular cusps. 

A special class of manifolds Let us denote by E g the closed orientable 
surface of genus g. For g ^ 2 and k ^ we define Ai g> k as the set of all 
compact orientable manifolds M having an ideal triangulation with g + k 
tetrahedra, and 

dM = £ g U f.L^J with Tj^ Ei. 

As a motivation for this definition, we mention here that an ideal triangu- 
lation of a manifold whose boundary is the union of E 9 and k tori contains 
at least g + k tetrahedra. So is the set of manifolds M having the 

smallest possible complexity, given the topological constraints on dM. 
Our first result shows that the class just introduced is very large: 

Proposition 5.1. • M. g ^k is non-empty precisely for g > k or g = k 

and g even; 

• The values of #M. g ,k f or small g and k are as shown in Table 1; 
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k = 


k = 1 


k = 2 


9 = 


2 


8 


1 


1 


9 = 


3 


74 


12 


1 


9 = 


4 


2340 


416 


51 


9 = 


5 


97568 


17900 


? 



Table 1: Some values of #M g ,k- 



• For any fixed k there exist constants C > c > such that 

g-9 < #Mg , k < g C 9 for g » 0. 

We now give the main statement we have about the elements of -M- g ,k and 
their topological and geometric invariants. We address the reader to [13] for 
the definition of the Heegaard genus of a triple (M, d M, d\M), and to [IE] 
for the definition of the Turaev-Viro invariants TV r (M) of a compact 3- 
manifold M, for r ^ 2. 

Theorem 5.2. Let M € M. g ^- The following holds: 

1. M\ (UTj) is hyperbolic, and its volume depends only on g and k; 

2. M has a unique ideal triangulation with g + k tetrahedra, which gives 
the canonical Kojima decomposition of M \ (UTj); 

3. M has complexity g + k; 

4- The Heegaard genus of (M, S g , UTj) is g + 1; 
5. #i(M;Z) = 

T/ie Turaev-Viro invariant TV r (M) o/M depends only on r, g and k. 

The first two points of this theorem are established using precisely the 
philosophy of moduli, equations, and tilts sketched in Section 0] Namely, 
given a minimal triangulation of M £ -M.g,ki ° ne shows that each tetrahedron 
has either one or no vertex asymptotic to a toric cusp. If there is one such 
vertex, one chooses the dihedral angles to be ir/3 at the edges ending at 
that vertex, and to be all equal to some a at the other edges. If there 
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is no such vertex, one chooses the tetrahedron to have all dihedral angles 
equal to some 0. Using a continuity argument then one sees that there exist 
values of a and f3 satisfying the consistency equations. Completeness and 
the computation of tilts is then straight-forward. 

We also notice that the previous results show the tremendous power 
of hyperbolic geometry compared to the topological invariants: the classes 
M g ,k are extremely large, and only hyperbolic geometry is able to distin- 
guish their elements from each other. 

We refrain from recalling the precise statements and details here, but we 
want to mention that a thorough analysis of the Dehn fillings of the elements 
of M g ,k was carried out in [Ej, leading in particular to the solution of a 
problem raised by Gordon and Wu [2U1 [551 ISH] on the maximal distance 
between a boundary-reducible and a non-acylindrical slope on a "large" 
hyperbolic manifold. 

Experimental results We describe here 7i n for n ^ 4. It is easy to see 
that is empty. Moreover, H.2 was shown in JHj to have 8 elements, all 
with the same volume ~ 6.451990. Therefore H.2 coincides with the set A^2,o 
discussed above, so all its members share the same invariants, except the 
hyperbolic structure itself. 

The strategy mentioned in Section [21 has been implemented in ^3] to 
classify H3 and Ti^, leading to the following results. We denote here by 
fC(Y) the Kojima canonical decomposition of a hyperbolic Y, and we recall 
that TL n is the set of candidate hyperbolic 3-manifolds of complexity n. We 
emphasize that the next results indeed have an experimental nature, but 
their validity was confirmed by a number of computations by hand and cross- 
checks. We also mention that all the values of the volume are approximate 
ones. More accurate approximations are available on the web |57| . 

Theorem 5.3. • H.3 coincides with H.3 and it has 151 elements; 

• -^3,0 C H.3 consists of 74 elements of volume 10.428602; 

• -M-2,1 C H.3 consists of a single manifold of volume 7.797637; 

• The 76 elements 0/H3 \ (A^o U A42,i) all have boundary T,2, and they 
split as follows: 

— 73 compact Y 's with IC(Y) consisting of three tetrahedra; the vol- 
ume function attains on them 15 different values ranging between 
7.107592 and 8.513926, with maximal multiplicity 9; 
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— 3 compact Y's with JC(Y) consisting of four tetrahedra; they all 
have the same volume 7.758268. 

Theorem 5.4. • Ti^ has 5033 elements, and and has 6 more; 

• Mifi C Ha has 2340 elements of volume 14.238170; 

• -^3,1 C Tii has 12 elements of volume 11.812681; 

• C Tii has a single element of volume 9.134475; 

• The 2680 elements Y of TCa not belonging to any hA g ^-g split as de- 
scribed in Table 2 according to the boundary ofY( columns ) and type 
of blocks of the Kojima decomposition (rows). Some volume informa- 
tion has also been inserted in each box, namely the unique value of 
volume if there is one, or the minimum and the maximum of volume, 
the number of values it attains, and the maximal multiplicity of these 
values. 

We conclude this section with some information on the actual computer 
implementation of the enumeration process. We solve the hyperbolicity 
equations using Newton's method with partial pivoting, after explicitly writ- 
ing the derivatives of the length function. Convergence to the solution is al- 
ways extremely fast, and it can be checked to be stable under modifications 
of the numerical parameters involved in the implementation of Newton's 
method. Concerning the Kojima decomposition, we mention that the evolu- 
tion of a triangulation toward the canonical decomposition is not quite sure 
to converge in general, but it always does in practice. We also point out that 
our computer program is only able to handle triangulations: whenever some 
mixed negative and zero tilts appear, the canonical decomposition must be 
worked out by hand and actually proved not to be a triangulation. 

6 Manifolds with annular cusps 

Given a hyperbolic manifold Y with annular cusps, one may decide to use 
either c(Y) of c(Y) as a definition of the complexity of Y, the former being 
more natural from a topological viewpoint, and the latter from a geometric 
viewpoint. However, as announced in Proposition 11.51 and explained in this 
section, neither definition allows to employ the powerful techniques of spe- 
cial spines and ideal triangulations to actually carry out the enumeration 
of manifolds in order of increasing complexity. For this reason the under- 
standing of manifolds with annular cusps is still very limited. This section 
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£3 


£ 2 


E2, 1 cusp 


4 tetrahedra 


1936 


555 


16 




min = 11.113262 


min = 7.378628 


min = 8.446655 




max = 12.903981 


max = 10.292422 


max = 9.774939 




values = 59 


values = 169 


values = 8 




max mult = 138 


max mult = 27 


max mult = 3 


5 tetrahedra 


42 

vol = 11.796442 


41 

min = 8.511458 
max = 9.719900 
values = 16 
max mult = 6 




6 tetrahedra 




3 

vol = 8.297977 




8 tetrahedra 




3 

vol = 8.572927 




1 octahedron 


56 


14 




(regular) 


vol = 11.448776 


vol = 9.415842 




1 octahedron 




8 




(non-regular) 




vol = 8.739252 




2 pyramids 




4 


2 


with square 




vol = 9.044841 


vol = 8.681738 


basis 









Table 2: Elements of H.4, not belonging to any A4 g A- g - 
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is only devoted to the description of a special class of such manifolds, but 
the properties of this class are already sufficient to show that the set of hy- 
perbolic manifolds with annular cusps is very large, and that there is very 
little control on the topology of the compactifications of such manifolds. We 
will also see that the geometric information on this special class allows to 
prove some very interesting (and apparently unrelated) results. We address 
the reader to [H] for all proofs and further details. 

Combinatorial triangulations We denote by T n the set of all possible 
simplicial pairings between the faces of n tetrahedra. We view two pairings 
to be the same if they coincide combinatorially, and we note that in general 
a pairing can fail give an ideal triangulation of a manifold, because the link 
of the midpoint of an edge can be the projective plane. Moreover, if a 
pairing actually gives a manifold, then this manifold need not be orientable. 
However, if we fix an orientation on the tetrahedron and require all the 
pairing maps to reverse the orientation, the result is an ideal triangulation 
of an orientable manifold. We denote by T the union of all 

Given T £ T we define a space Y(T) by performing the pairings in T 
and then removing first an open regular neighbourhood of the vertices and 
then a closed regular neighbourhood of the edges. 

Remark 6.1. The manifold Y(T) is a non-compact one with boundary, and 
its topological ends have the shape of the product of an annulus and a closed 
half-line. Let us define Y(T) and Y(T) as the natural compactifications of 
Y(T), as we did in Section ^ for a hyperbolic Y. Then Y(T) is a (possibly 
non-orientable) handlebody, with genus n + 1 if T G T n . Moreover, if T is 
an ideal triangulation of a manifold M, then Y(T) is homeomorphic to M. 

Relative handlebodies and hyperbolicity For T € T the manifold 
Y(T) can be viewed as H \ T, where H is a handlebody and T C dH is a 
system of disjoint loops. At the risk of a little ambiguity, which eventually 
will not create any serious problem, we then identify Y(T) to the pair (H, T), 
and, following Johannson 23= , we call such a pair a relative handlebody. We 
denote by A n the set of all y(T)'s as T varies in T n , and by A the union of 
all A n 's. 

To state our first result we call complexity of a relative handlebody (H, T) 
the minimum of |rn<9-D| where D is a system of disjoint properly embedded 
discs in H such that dD cuts dH into a union of pairs of pants. 
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Proposition 6.2. • For allT G T £/te relative handlebody Y(T) is hy- 
perbolic. 

• The map T i— > V(T) gives a bijection between T and A. 

• Among the hyperbolic (H, T) 's of genus n+1, the elements of An can be 
characterized as those having minimal complexity, equal to 10 • n, and 
as those having minimal volume, equal to n ■ vo, where vo ~ 3.66386 
is the volume of a hyperbolic regular ideal octahedron. 

The proof of this result entirely depends on the techniques described 
in Section |SJ One first notices that a regular ultra-ideal tetrahedron with 
all dihedral angles equal to is actually a regular ideal octahedron with a 
checkerboard coloring of the faces, the white ones being internal faces and 
the black ones being truncation faces. Then one sees that for each given 
T £ T n the corresponding Y(T) can be constructed by gluing according to 
T the white faces of n regular ideal octahedra. To conclude, one proves that 
T itself is the canonical Kojima decomposition of Y{T). 

Remark 16.11 and Proposition 16.21 are already sufficient to prove Propo- 
sition 11.51 which shows, as already mentioned, that there are very many 
manifolds with annular cusps and that their topology is rather arbitrary. 
The next result goes in the same direction. We call tangle in a compact 
3-manifold with boundary a finite union of disjoint properly embedded arcs. 

Corollary 6.3. Every tangle in every compact 3-manifold is contained in a 
tangle whose complement lies in A n for some n (in particular, it is hyperbolic 
with annular cusps). 

An unexpected consequence of Proposition 16.21 is also the following: 

Theorem 6.4. Let Tq and T\ be triangulations of the same compact 3- 
manifold, with finite and/or ideal vertices, possibly with multiple and self- 
adjacencies. Assume that the 1-skeleta o/To and Ti coincide. Then Tq and 
T\ are isotopic relatively to the 1-skeleton. 

Doubles and Dehn filling We state in this paragraph some of the sur- 
prising results one can establish starting from the construction of the hy- 
perbolic relative handlebodies Y{T). 

For T E T, we define D(T) to be the "orientable double" of Y(T), namely 
either the union of two copies of Y(T) along the boundary, when Y(T) is 
orientable, or the quotient of the orientation covering of Y(T) under the 
restriction to the boundary of the involution, when Y{T) is non-orient able. 



23 



We define T> as the set of all Y(T) for T € T. The next result shows that the 
family T> is completely classified in combinatorial terms and it is universal 
for 3-manifolds under the operation of Dehn filling: 

Theorem 6.5. • Every member ofT> is an orientable cusped hyperbolic 
manifold without boundary, and it is the complement of a link in a 
connected sum of some copies of S 2 x S ; 

• The correspondence T i— > D(T) defines a bisection between T and T>. 

• Every closed orientable 3-manifold is a Dehn filling of a manifold in 
V. 

Since a Dehn filling of a hyperbolic manifold "typically" is hyperbolic, 
the class T> provides a powerful method to construct closed hyperbolic man- 
ifolds. As an application of this method one can establish the following 
refinement of the main result of |25| : 

Proposition 6.6. There exists c > such that, given a finite group G, 
there is a closed orientable hyperbolic 3-manifold M with Isom(M) = G and 
Vol(M) < c- |Gf. 

We now recall that, according to Thurston's hyperbolic Dehn filling the- 
orem, on each cusp of a finite-volume hyperbolic 3-manifold there is only 
a finite number of slopes filling along which one gets a non-hyperbolic 3- 
manifold. These slopes are called exceptional, and a considerable effort has 
been devoted to understanding them . If T is a triangulation, the hy- 
perbolic manifold D(T) has a preferred horospherical cusp section, and each 
component of this section corresponds to an edge of T. Moreover the valence 
of the edge gives a lower bound for the length of the second shortest geodesic 
on the component. This fact and the Agol-Lackenby 6-theorem pflEE] imply 
the following: 

Proposition 6.7. If every edge of T has valence at least 7 then there is at 
most one exceptional slope on each cusp of D(T). 

We conclude this paper by showing that in some cases the combinatorics 
of a triangulation is already sufficient to prove hyperbolicity of the under- 
lying manifolds. To motivate our result, suppose first that a 3-manifold M 
has an ideal triangulation T in which each edge has valence at least 6. An 
easy argument shows that x(T) ^ 0, and that x{T) = precisely when all 
valences are 6. Moreover, in the last case, the boundary of M is a disjoint 
union of tori and Klein bottles, and Thurston's hyperbolicity equations for 



24 



cusped manifolds have a very simple solution, given by regular ideal tetra- 
hedra. Analogously, if all edges have one and the same valence v ^ 7, the 
hyperbolicity equations for the geodesic boundary case have a simple solu- 
tion, given by regular truncated tetrahedra with dihedral angles 2tt/v. An 
argument based on the Agol-Lackenby machinery ^ I28| allows to generalize 
these facts as follows: 

Proposition 6.8. // M has an ideal triangulation T whose edges have va- 
lence at least 6, then M is hyperbolic, and the edges of T are homotopically 
non-trivial relative to dM . 

According to the last assertion of this result, the edges of T can be 
straightened to geodesies in M without entirely disappearing into infinity, 
which suggests that T itself can be straightened. We address the reader 
to |41| for more information on the general problem of existence of a trian- 
gulation which can be straightened to an ideal one. 
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